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Abstract 

Background: Influenza viruses can modulate and hijack several cellular signalling pathways to efficiently support 
their replication. We recently investigated and compared the cellular gene expression profiles of human lung A549 
cells infected by five different subtypes of human and avian influenza viruses (Josset et al. Plos One 2010). Using 
these transcriptomic data, we have focused our analysis on the modulation of the p53 pathway in response to 
influenza infection. 

Results: Our results were supported by both RT-qPCR and western blot analyses and reveal multiple alterations of 
the p53 pathway during infection. A down-regulation of mRNA expression was observed for the main regulators of 
p53 protein stability during infection by the complete set of viruses tested, and a significant decrease in p53 mRNA 
expression was also observed in H5N1 infected cells. In addition, several p53 target genes were also down- 
regulated by these influenza viruses and the expression of their product reduced. 

Conclusions: Our data reveal that influenza viruses cause an overall down-regulation of the host p53 pathway and 
highlight this pathway and p53 protein itself as important viral targets in the altering of apoptotic processes and in 
cell-cycle regulation. 



Background 

Influenza viruses belong to the Orthomyxoviridae family 
of enveloped viruses containing a segmented genome of 
single stranded negative RNA. Among the three influ- 
enza types (A, B and C), type A is the most virulent 
pathogen with a diversity represented by the combina- 
tion of 16 H and 9 N different subtypes (e.g. H1N1, 
H5N1) [1]. Influenza A viruses are the most serious 
threat to public health, with the potential to cause glo- 
bal pandemics as was illustrated in 2009 with the emer- 
gence of H1N1 SOIV [2]. 

All known subtypes of the influenza A virus are main- 
tained in wild waterfowl, the natural reservoir of these 
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viruses [3]. Current human circulating influenza A 
subtypes are H1N1 and H3N2. While extensive viral 
diversity is responsible for the subtype-specific virus-host 
interactions, many common functional features are also 
shared among viruses. Influenza infection alters host cel- 
lular homeostasis via the combination of the virally- 
induced alteration of biological machineries/pathways and 
the cellular antiviral response triggered by intracellular 
signalling cascades. Influenza viruses are able to activate/ 
inhibit and hijack several cellular signalling pathways to 
efficiently support their own replication [4]. 

The development of high-throughput 'omic' studies 
has increased our understanding of viral-host interac- 
tions. Numerous studies have described host gene 
expression modifications induced upon viral infection 
in vitro, in animal models and in patients [5-9]. For exam- 
ple, a global mRNA profiling study of MDCK-infected 
cells has revealed an important role of the NF-kappaB 
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signalling pathway for the H5N1 subtype and to a lesser 
extent for the currently circulating H3N2 strain [10]. 
Other intracellular signalling cascades are also induced by 
infection, in particular the mitogen-activated protein 
kinase (MAPK) and the PI3K/Akt pathways, both of which 
activate downstream transcription factors thereby affecting 
host gene expression [11,12]. While a small number of 
transcriptional studies have compared the cellular 
response induced by infection with highly pathogenic 
strains such as H1N1 1918 and H5N1 to that induced 
with low pathogenic strains, a systematic analysis and 
comparison of host cell mRNA expression during infec- 
tion by several genetically diverging influenza subtypes has 
not yet been performed [7,8,13,14]. 

Using a nylon microarray, we recently characterized 
common gene expression changes induced during the 
infection of human lung epithelial cells by five different 
influenza A viruses. This shared signature was exploited 
to find new molecules that act on host metabolic path- 
ways to bring about an antiviral effect against several 
subtypes [15]. Moreover, our systematic transcriptomic 
study also provided a gene expression database with 
which we compared the cellular responses induced by 
different influenza viruses. In the present study, we used 
the same data set to analyze strain specificity in regard 
to several host pathways. We observed that each viral 
signature was mainly associated with cellular pathways 
linked to cell death and cell cycle progression (i.e. cell 
growth). At the crossroads of all of these cellular pro- 
cesses, the p53 pathway became the focus of our investi- 
gation. Our analyses, validated by both RT-qPCR and 
western blot analyses, reveal multiple alterations within 
the p53 pathway during influenza A infection. 

Methods 

Cell lines, viruses and infection 

Human lung epithelial A549 cells (ATCC CCL-185) or 
human colon carcinoma HCT116 cell lines were grown as 
monolayers in Dulbecco's modified Eagle's medium 
(DMEM, Gibco), supplemented with 10% heat-inactivated 
foetal bovine serum, 2 mM L-glutamine, 100 U/mL of 
penicillin and 100 mg/mL of streptomycin sulphate, at 
37°C. The HCT116 cell lines containing a p53 wild-type 
(p53+/+) and a p53-deleted derivative (p53-/-) were 
gifts from Dr. Bert Vogelstein (John Hopkins University, 
Baltimore, MD, USA) [16]. Influenza viruses A/New Cale- 
donia/20/99 (H1N1), A/Moscow/10/99 (H3N2), A/Lyon/ 
969/09 (H1N1 SOIV), A/Turkey/582/2006 (H5N1), 
A/Finch/England/2051/94 (H5N2), and A/Chicken/Italy/ 
2076/99 (H7N1) were produced in MDCK cells in EMEM 
supplemented with 2 mM L-glutamine, 100 U/mL of peni- 
cillin, 100 mg/mL of streptomycin sulphate and 1 mg/mL 
of trypsin. Viruses were titrated to determine the 50% 
tissue culture infective dose (TCID50) in MDCK cells as 



described in our previous study [17]. Confluent A549 cells 
were infected with influenza viruses at a MOI of 1 or 10-3 
for one hour in a minimal volume of DMEM supplemen- 
ted with 2 mM L-glutamine, 100 U of penicillin/mL, 
100 ug of streptomycin sulphate/mL and 0.5 ug of trypsin/ 
mL (infection medium) at 37°C. The cells were then over- 
laid with fresh infection medium and incubated at 37°C 
for 24 h. Viral kinetics for the six different influenza 
viruses have been determined previously [15]. No markers 
of apoptosis were detected at 24h post infection (MOI 1) 
by western blot (capspase 3 cleavage) and no signs of cell 
death were observed (data not shown). 

Microarray analysis 

Gene expression data of A549 cells infected by H1N1, 
H3N2, H5N1, H5N2 or H7N1 influenza viruses or 
mock infected (five replicates in each group) were 
obtained in our previous study [15]. Data sets are avail- 
able publicly from the Gene Expression Omnibus (GEO) 
database (http://www.ncbi.nlm.nih.gov/geo/) under 
accession number GSE22319. Five supervised analyses 
between groups of infected versus mock samples were 
conducted using the Significance Analysis of Microarray 
algorithm (SAM) [18], and the siggenes library (vl.18.0) 
[19]. Ingenuity Pathway Analysis 5.0 (IPA) (Ingenuity 
Systems, Redwood City, CA, USA) was used to select, 
annotate, and visualize genes according to function and 
pathway (Gene Ontology). Additional gene annotation 
was provided by the Interferome Database [20] and 
DAVID database (http://david.abcc.ncifcrf.gov/). Heat 
maps were produced by the heat map function from R 
that uses hierarchical clustering with Euclidean distance 
metric and the complete linkage method to generate the 
hierarchical tree (http://www.r-project.org/). 

Validation by RT-qPCR 

Total RNAs were extracted using the RNAeasy Mini Kit 
(Qiagen). Reverse-transcription was performed on 1 ug 
of total RNAs using the Superscript II enzyme (Invitro- 
gen) at 42°C. Quantification of the level of different 
mRNAs of interest was performed by real-time PCR on 
a MX3005P apparatus (Stratagene). Briefly, 20 ng of 
cDNAs were amplified using 0.8 uM of each primer, 
0.4 uM of probe (cf. Table 1) and IX Taqman Universal 
Master Mix (Applied Biosystems). All data were normal- 
ized to the internal standard Actin. For each single-well 
amplification reaction, a threshold cycle (Ct) was 
observed in the exponential phase of amplification. Rela- 
tive changes in gene expression were determined using 
the 2AACt method and reported as the n-fold difference 
relative to a control cDNA (mock, uninfected cells) pre- 
pared in parallel with the experimental cDNAs (infected 
cells) [21]. The mRNA levels were measured in triplicate 
in two independent experiments. Quantification of M 
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Table 1 List of primers used in this study 



Primers & probes 


Sequence 5'-3' 


p53 (Exon8-9) 


Forward 


f~ A A f~ A f~ A AT" f~f~f~ Aflf AAA f~ f~ 

LAA LAL AA I LIC LLL AAL AAA LL 




Reverse 


I LL AIL LAG I LL I I I L I I L I I I L 




Probe 


a /~r^ a f — r a a f f~r~ a rr a f — vc e~t~r~ a at 
ALL AL I AAL LLA LLA L I L LLL AAL 


p21 


Forward 


r~ a t~T( — rc a (~ct — rcr~ a a a a ccc 
LAL I L I LALLL I LLAAAALLL 




Reverse 


LLLLA I I ALLLL I I LL I L I I 




Probe 


(~~T A f~f~ A ( — Yf~f~ AAA C f~ Cf~f~ f~ t — Tf~ AT/ — T 

L I ALLAL I LLAAALLLLLLL I LA I L I 


Rav 
DdX 


Forward 


AL I LLLLLLLALALL I L I I 




Reverse 


CC AAAfTTAfiAAA ACCCnCCAC AA 




Probe 


GAGC I GALA I G I I I I L I GACGGLAAL I I CAAL I 


Bcl-XL 


Forward 


TCG TTG TCT AGG CTT TCC ACG 




Reverse 


GGT CGG ATT CTC GGG TTT 




Probe 


ACA GTG CCC CGC CGA AGG AGA 


Influenza M 


Forward 


CTTCTAACCGAGGTCGAAACGTA 




Reverse 


GGTGACAGGATTGGTCTTGTCTTTA 




Probe 


TCAGGCCCCCTCAAAGCCGAG 



viral genome copies released in supernatants was per- 
formed according to previously published work [22]. 

Western blot analysis 

Total proteins were extracted by scraping and syringing 
cells in IX NuPAGE LDS buffer (Invitrogen). Fifteen to 
thirteen micrograms of total proteins were then sepa- 
rated on pre-cast 10% NuPAGE gels (Invitrogen). To 
detect the different proteins of interest, the following 
antibodies were used: mouse monoclonal DO-1 (p53) 
and SX118 (p21, Cell signalling), Rabbit polyclonals Bax, 
Akt and Bcl-XL (Cell signaling). The mouse monoclonal 
Ku80 antibody was used as a loading control (Abeam). 
The relative protein levels (RPL) were calculated by den- 
sitometry analysis performed with the help of the ImageJ 
software (http://rsbweb.nih.gov/ij/). 

Results 

Specific viral signatures and commonly targeted host 
genes 

We have previously reported variations in gene expres- 
sion at 24 h post-infection in human pulmonary epithe- 
lial A549 cells infected with different subtypes of 
influenza viruses: human A/New Caledonia/20/99 
(H1N1), A/Moscow/10/99 (H3N2), avian A/Turkey/582/ 
2006 (H5N1), A/Finch/England/2051/94 (H5N2), A/ 
Chicken/Italy/2076/99 (H7N1) influenza viral strains or 
mock control infection (five replicates in each group) 
[15]. Compared to the mock control infection, changes 
were observed in 36 cellular genes for H1N1, 2298 
genes for H3N2, 1455 genes for H5N1, 1510 genes for 
H5N2 and 3020 genes for H7N1 (SAM algorithm with a 
False Discovery Rate of 10%). The small changes 
induced upon H1N1 infection were correlated to its low 



level of replication in our experimental conditions [15]. 
While 83% of the host genes altered following H1N1 
infection were over-expressed, a similar number of 
genes were over- and under-expressed for each of the 
other viruses (figure 1A). Furthermore, we found that 
the expression levels for a total of 300 genes (represent- 
ing 3.4% of the genes considered present on the chip) 
differed significantly between mock and all infected 
samples including H1N1 [15]. 

To determine the specific cellular genes showing 
altered expression within each of the different strain- 
specific signatures, we built the Venn diagram shown in 
figure IB. Common changes in the expression of 
numerous genes (1007) were induced by three of the 
viruses, H3N2, H5N2 and H7N1. The expression of 705 
of these genes was also altered when four viruses were 
considered: H3N2, H5N1, H5N2 and H7N1. Since 
H1N1 deregulated only a small number of genes, the 
number of those the expression of which was modified 
by all strains was limited to 20. Their associated biologi- 
cal responses were described in our previous study and 
are mainly linked to the immune response in accor- 
dance with other investigations [15]. 

Biological processes associated with each viral signature 

We then determined the Gene Ontology (GO) biological 
processes (or terms) significantly associated with each 
viral signature using the DAVID database. With a 
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H1N1 


H3N2 HSN1 


H5N2 


H7N1 


Total 
UP 
DOWN 
IRGs 


36 

30 (83,3%) 
6(16,7%) 
17(46%) 


2296 1455 
1098(47.8%) 660(45,4%) 
1200(52,2%) 794(54.6%) 
351(15.3%) 256(17.6%) 


1510 3020 
687(45,5%) 1541(51.1%) 
623(54,5%) 1479(48,9%) 
231(15.3%) 456(15,1%) 
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Figure 1 Summary of genes expressed in response to H1N1, 
H3N2, H5N1, H5N2 and H7N1 viral infection A Genes 
significantly regulated (SAM, FDR = 10%) in response to the 
different influenza viruses compared to mock infection at 24 hpi are 
shown. B. Venn-diagram showing the genes co-regulated by several 
viruses 
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p-value < 0.01, no one term was shared among the five 
viruses. The 36 genes deregulated upon H1N1 infection 
encode proteins involved in the cellular immune 
response, transcription, biogenesis of protein complexes 
and nucleic acid metabolic processes. In contrast, the 
signatures of H3N2, H5N1, H5N2 and H7N1 viruses 
were each associated with 33 common terms (p-value 
< 0.01) (figure 2), covering all the biological processes of 
DAVID database. Two major classes of biological pro- 
cesses could be distinguished. Firstly, there are those 
that are affected by the four viral strains through the 
modulation of a common gene subset. For example, 
each of the viral strains modulated the expression of 10 
common genes involved in the glucose catabolic pro- 
cess", while only a few genes involved in this process 
were differentially modulated in a strain-dependent 
manner. Secondly, there are those that are affected by 
each of the four viral strains through the modulation of 
distinct gene subsets in a strain-dependent manner. For 
example, while 35 genes all associated with the "apopto- 
sis" process were modulated by the four viruses, 11 of 
these were specifically modified upon H3N2 infection, 
11 by H5N1, 7 by H5N2 and 35 by H7N1. 

Among the 33 common terms associated with the 
viral signatures, more than 60% (n = 20) were related to 
the control of cell death or cell cycle progression. Other 
common biological processes were related to angiogen- 
esis, protein complex biogenesis or metabolism [15]. 
Interestingly, despite the abovementioned current asso- 
ciation made between influenza virus-induced host cell 
responses and the cellular immune response, we found 



no terms shared by all viruses relating to the immune 
response. Moreover, by consulting the whole list of bio- 
logical processes associated to each virus signature we 
found the term "immune response" significantly mapped 
to each viral profile, except for the H5N2 virus. How- 
ever, the comparison of our differentially expressed gene 
lists to Interferome, an IFN-regulated gene database, 
revealed that a similar proportion of genes regulated 
upon infection by H3N2, H5N1, H5N2 or H7N1 
(around 15%) was related to the IFN response, whereas 
48% of the genes responding to H1N1 infection were 
regulated by IFN [20]. 

In summary, analysis of our data showed that more 
than half of the genes responding to infection by H3N2, 
H5N1, H5N2 and H7N1 viruses, are involved in cell 
death, cell cycle progression or the cellular immune 
response. All these biological processes are regulated by 
a central signalling pathway, the p53 pathway [23]. 
Indeed, in addition to regulating the expression of genes 
involved in apoptosis and cell cycle arrest, p53 also reg- 
ulates the expression of type I IFN as well as several 
genes carrying IFN-stimulated response elements [24]. 
For these reasons, and in light of several proposals of 
p53 being a functional interactant of influenza [25-27], 
we focused our investigations on modifications of the 
p53 pathway in response to influenza infection. 

Infection induces a transcriptional down-regulation of the 
upstream signal part of the p53 pathway 

In response to stress, such as viral infection, several sig- 
nalling mediators relay the stress stimuli to the p53 





Total Numbei of genes In each vital signature (p value) 


common 


Number ol 


qenes 


unique to each sign 


jtuie 


Tmr 


H3N2 


H4N1 


H5N2 


H7N1 


gone* 


K3N2 


ir.Ni 


H5H2 


H7N1 


GO 0001558-toguialion of coll growth 

GO 0001568-blood vessel development 


51 (B 7o-05) 

67 (1 40-06) 


38(5 0e-05) 
43 (2 Oe-04) 


33 (BOe-05) 

34(2 5*03) 


63 P 80-05) 

85(6 10-08) 


18 

12 


5 
11 


5 
1 




12 
18 


GO 0001944-vasculalure development 


68 (1 60-06) 


43 (3 5*04) 


34 p 7*03) 


87 (4 4e-08) 


t: 


11 


1 




13 


GO 0006007-glucose catabolic process 


19 (1 Be-03) 


13 (9 1e-03) 


13(2 0e-03) 


22 (2 7*03) 


to 


2 


0 




4 


GO 0006461 -protein complei assembly 


104 (8 90-04) 


72 (9 7e-04) 


61 (1 6e-03) 


139(2 0*05) 


27 


13 


7 




33 


GO 0006915- apoptosis 


124 (2 8e-04) 


83(1 le-03) 


69(2 6003) 


158 (8 1*05) 


35 


11 


11 




36 


GO 0006916- anti- apoptosis 


57(6 30-06) 


32 (9 70-03) 


32 (5 50-04) 


68 (1 Oe-05) 


15 


6 


4 




10 


GO00O6917~induction of apoptosis 


77(3 40-05) 


50 (9 5e-04) 


42(21o03) 


10E (2 0*08) 


18 


6 


7 




26 


GO 0007010-cytoskeleton organization 


93 (5 30-04) 


77 p 3*07) 


54 (1 7*03) 


1 18 (1 9*04) 


30 


8 


10 




21 


GO 0007049-cell cycle 

GO 0007243-protem kinase cascade 


168(9 90-07) 

83 (2 le-04) 


125 (1 3e-03) 

55 (1 6.-03) 


84(4 7e-03) 

60p0e-07) 


206(2 6*06) 
117(6 5*08) 


41 

26 


21 
9 


17 

5 




39 
19 


GO 0008219-cell death 


149(4 6 0-05) 


101 (t 5*04) 


81 (1 9*03) 


192(4 0*06) 


41 


14 


12 




45 


GO 0008361 -regulation of celt size 


53 (1 10-04) 


36 (7 9e-04) 


32 (5 5e-04) 


61 (81*04) 


18 


6 


5 




11 


GO 0010941 -regulation of cell death 


202 (2 7o-13) 


128(4 1*08) 


113(1 1*08) 


255 (1 3*15) 


57 


18 


18 




44 


GO 001D942-posrtive regulation of cell death 


111 (2 Oe-08) 


71 (1 7e-05) 


59 p 9e-05) 


147(5 2*12) 


31 


9 


8 




33 


GO 0012501 -programmed cell death 


127 (1 6e-04) 


86 (4 5e-04) 


71 (1 7*03) 


162 p 7*05) 


37 


11 


11 




37 


GO 0012502-induction of programmed cell death 


77(3 8045) 


50(1 00-03) 


42(2 3*03) 


106 (2 5*08) 


18 


6 


7 




26 


GO 0016265-death 


150 (4 5e-05) 


104 (4 7*05) 


82(1 5e-03) 


193 (4 3*06) 


42 


14 


14 




45 


GO 0030155-regulation of cell adhesion 


35(2 20-03) 


27 tj 8e-04) 


24 (6 Oe-04) 


42 P Oe-03) 


13 


7 


5 




8 


GO 0032535-regulation of cellular component size 


65(1 6*04) 


40 (8 6e-03) 


37(2 1*03) 


78 p 8*04) 


16 


7 


S 




IE 


GO 0033554-cellular response to stress 


118(2 4e-04) 


85 (5 Se-05) 


72 (1 2*04) 


150(7 6*06) 


38 


14 


15 




32 


GO 0040006-regulation of growth 


88 (4 1e-07) 


58 p 4*05) 


53(5 6*06) 


101 (1 4*05) 


27 


10 


10 




18 


GO 0042361 -regulation of apoptosis 


138 p 3*13) 


127 p 3*08) 


112 p 4*09) 


250(59*15) 


56 


18 


16 




43 


GO 0043065-posrtive regulation of apoptosis 


110 (2 1e-08) 


71 (1 1*05) 


59(7 1*05) 


146 (4 1*12) 


31 


9 


8 




33 


GO 0043066-negatrve regulation of apoptosis 


86 (7 7e-06) 


62 (2 9*03) 


51 (6 9*05) 


109 (8 2*07) 


20 


8 


10 




15 


GO 0043067-regulation of programmed cell death 


201 (3 6o-13) 


128 p 3*08) 


113(8 5*09) 


253 P 0*15) 


57 


18 


18 




43 


GO 0043068-positrve regulation of programmed cell death 


110 (3 10-08) 


71 (1 4e-05) 


59 (8 7*05) 


146(7 4*12) 


31 


9 


8 




33 


GO 0043069-negatrve regulation of programmed cell death 


88 (4 le-06) 


52 p 9o03) 


51 (1 0*04) 


111 (51*07) 


20 


8 


10 


S 


IS 


GO 0043933-macromolecular complex subunrt organization 


146 (8 2e-0S) 


92 p 7e-03) 


76 0 4*03) 


187 (1 4*05) 


35 


20 


10 


5 


45 


GO 0046907-mtracellular transport 


129 (1 5e-03) 


84 (8 2e-03) 


73 (4 6*03) 


171 (6 5*05) 


X 


12 


3 


5 


■;: 


GO 0048514-blood vessel morphogenesis 


59 p Oe-06) 


39 (1 4*04) 


30 P 3*03) 


77(2 5*08) 


11 


9 


1 


3 


17 


GO 0060548-negatrve regulation of cell death 


88 (4 6e-06) 


52(4 le-03) 


51 (1 1*04) 


112 p 2*07) 


20 


8 


10 


S 


16 


GO 0070271 -protein complex biogenesis 


104 (8 9e-04) 


72 p 7*04) 


61 (1 6*03) 


139 P 0*05) 


27 


13 


7 


5 


33 



Figure 2 Significantly enriched GO Biological Processes (p-value < 0.01) common to H3N2, H5N1, H5N2 and H7N1 specific signatures 

Terms related to the p53 pathway are highlighted in grey. 
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protein, which is consequently stabilized and activated 
mainly by post-translational modifications including 
phosphorylations and acetylations [28]. Such modifica- 
tions disrupt interactions between p53 and the ubiquitin 
ligase Hdm2, responsible for the rapid degradation of 
p53 in basal conditions [29,30]. Once activated, p53 reg- 
ulates cell cycle arrest and apoptosis through the direct 
modulation of host cellular gene expression. 

We firstly focused our attention on the upstream sig- 
nal part of the p53 pathway. According to our data, all 
the main kinases known to directly regulate p53 phos- 
phorylation status (DNA-PK, Chk2, JNK1 and Gsk3|3) 
are under-expressed following infection by H3N2, 
H5N1, H5N2 and H7N1 viruses (figure 3). For example, 
n-fold DNA-PK mRNA changes of -1.78, -2.12, -1.41 
and -1.56 were measured for H3N2, H5N1, H5N2 and 
H7N1, respectively (figure 3). We also focused our 
attention on the phosphatidylinositol 3-kinase (PI3K)/ 
Akt pathway, which directly regulates Hdm2 expression 



and thus endogenous p53 activity [31]. Interestingly, our 
results indicate a decreased expression of both AKT1 
and its negative regulator PTEN (figure 3). These results 
are in accordance with the activation of the phosphatidy- 
linositol 3-kinase (PI3K)/Akt pathway during influenza 
infection, as already described in the literature [32] . 

Altogether, our results reveal that influenza A infection 
induces the transcriptional down-regulation of several 
genes encoding host factors belonging to the upstream 
signal part of the p53 pathway and which act on the 
endogenous activity of p53. 

Differential effects of influenza viruses on p53 mRNA and 
protein levels 

Secondly, we assessed whether viral infection could also 
directly affect p53 mRNA expression. Interestingly, our 
data indicates that only one of the investigated viruses, 
H5N1, led to a significant decrease in the level of p53 
mRNA expression, with a fold change of -1.54 (figure 3). 




Figure 3 Transcriptional regulation of the genes belonging to the p53 signalling pathway (IPA canonical pathway) during influenza 
viral infection. Mean fold changes of the genes encoding regulators and targets of p53 that were differentially expressed by at least one virus, 
calculated as Lx^fMean expression levels in infected samples) - Log 2 (Mean expression levels in mock samples) and overlaid on the pathway in 
green for downregulated genes (FC < 0), and red for upregulated genes (FC > 0). Fold changes in each group of samples are depicted in the 
heat maps and were used to cluster hierarchically both samples and genes. Genes considered significant for each virus with the SAM procedure 
are indicated by a white X on the heat maps. 
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To assess such differential data, we infected another set 
of A549 cells with the different viruses used in our initial 
microarray study, including the low replicative H1N1 
virus to which we added the recent pandemic H1N1 
SOIV as an efficient productive H1N1 virus [28]. Total 
cell RNA was extracted 24 hours post-infection, similarly 
to that performed for the microarray analysis. The quan- 
tification of p53 mRNA levels was carried out by RT- 
qPCR using primers matching several different EST 
sequences used in the microarray chip. In accordance 
with our previous data, no significant change in the p53 
mRNA ratio was detected with H3N2, H5N2 and H7N1, 
nor with H1N1 or H1N1 SOIV, while a significant ratio 
change of 0.70 was observed in H5N1 infected cells (p < 
0.0005, figure 4A). The calculated fold change was -1.43, 
indicating the specific down-regulation of p53 mRNA 
expression in the context of H5N1 infection. 

We further investigated whether such influenza- 
induced variations in p53 mRNA levels could be corre- 
lated or not with changes in the endogenous p53 
protein level. For that, we reproduced the same experi- 
mental conditions as described above and performed 
western blots on infected cellular extracts. Compared to 
mock-infected cells, the p53 protein levels (Relative 
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Figure 4 A. Validation of microarray data by RT-qPCR The 

measured ratio change of p53 mRNA levels was subject to statistical 
analysis (student t-test, *** p-value < 0.001). B. p53 protein levels 
revealed by western blot. C. Relative protein levels (RPL) measured 
by densitometry analysis. 



Protein Level, RPL) measured by densitometry analysis 
were significantly increased in cells infected by all the 
viruses studied, with a 1.78 to 2.72 fold increase for 
respectively H7N1 and H1N1 (SOIV), and a smaller 
1.22 fold increase for H5N1 (figure 4B and 4C). 

Our data further support the previously reported non- 
correlation between p53 mRNA and protein levels in 
response to influenza infection, and are also in line with 
several other investigations reporting such a mismatch 
in the context of other cellular stresses [33]. Moreover, 
our results reveal an interesting contrast between H5N1 
and the other viruses. Indeed, H5N1 infection directly 
affects the p53 mRNA expression level, accounting for 
the attenuated increase in the p53 protein level observed 
in the H5N1 infected cell. 

Transcriptional down-regulation of downstream signalling 
of the p53 pathway upon infection 

Having determined an increase in p53 protein levels 
during the time-course of infection, we then focused 
our attention on the expression of p53 downstream tar- 
get genes. While several studies have previously shown 
an increase in p53 transcriptional activity during influ- 
enza virus infection by luciferase assay or the detection 
of phosphorylated p53 in infected cells [25,34], we 
observed in our experimental conditions a significant 
decrease in the transcriptional expression of several 
p53-targets genes during infection by H3N2, H5N1, 
H5N2 and H7N1 viruses (figure 3). These down-regu- 
lated genes include p21, 14-3-3, PERP, FAS, DR4/5, 
PIG3, BAX, Bcl-XL, PAI-1 and PCNA. For p21 
(CDKN1A) for example, the n-fold mRNA changes mea- 
sured were -1.81, -2.0 and -1.85 for H3N2, H5N1 and 
H7N1 infections, respectively. The expression of only 
two p53-target genes investigated, Teap (TP53INP1), 
and GADD45G was increased upon infection by the 4 
viruses (figure 3). The GADD45G n-fold mRNA changes 
measured were 1.61, 1.70, 1.26 and 1.60 for H3N2, 
H5N1, H5N2 and H7N1, respectively. It should be 
noted that TP53INP1 is also known to be regulated by 
p73, a p53 homologous protein, and must therefore be 
also regulated in a p53-independent manner [35]. 

To support these results, another set of A549 cells 
were infected with the broad set of viruses for 24 hpi 
and variations of endogenous p21, Bax and Bcl-XL 
mRNA expressions and protein levels were analyzed in 
parallel by RT-qPCR and western blot, respectively 
(figure 5 A, B and 5C). For p21, the mRNA level was 
significantly reduced only in H5Nl-infected samples 
(p < 0.05, mRNA ratio around 0.2), while the protein 
levels were reduced in cells infected by almost all the 
viruses studied (figure 5A). Moreover, a down-regulation 
of Bax expression was also clearly observed with signifi- 
cantly reduced mRNA levels in all infected cells, despite 
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Figure 5 Validation of microarray data by RT-qPCR and western blot: p53 target genes, (A) p21, (B) Bax, (C)Bcl-XL. 

correspond to p-values <0.05, <0,005, <0.001, respectively). 



and 



only minor changes at the protein level (figure 5B). We BAX and Bcl-XL, was considerably greater in the 



also confirmed the down-regulation of Bcl-XL expres- 
sion in cells infected by all viruses except H1N1 SOIV 
and H5N2 (figure 5C). At the protein level, Bcl-XL pro- 
tein level was decreased by most viruses. Of note, the 
down-regulation of the three p53-target genes, p21, 



context of H5N1 infection than with the other influenza 
A viruses. 

Altogether, our results reveal that influenza A infec- 
tion leads to the down-regulation of several p53 tran- 
scriptional targets, at both the mRNA and protein levels 
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for some, despite the relative increase in p53 protein 
levels common for all viruses. 

p53 contributes to an antiviral cellular state 

We then wished to assess whether such a global down- 
regulation of the p53 pathway could be associated with 
a functional effect with regards to viral replication. For 
this, we compared the production yield of one represen- 
tative influenza virus (subtype H3N2, circulating strain 
A/Moscow/ 10/99) in the HCT116 isogenic cell models 
(figure 6). HCT116 p53+/+ and p53-/- cells were 
infected with H3N2 at a MOI of 10-3 and the released 
viral genome copies were extracted from culture super- 
natants and quantified by RTqPCR at 24 and 48 hpi, as 
described in materials and methods. The number of 
genome copies reflecting viral production was 10 and 
100 fold higher in HCT116 p53-/- than in HCT116 p53 
+ /+ cells at 24 and 48h post-infection, respectively. 
These results show the effectiveness of p53 protein at 
slowing down viral production. In agreement with pre- 
vious reports [6,25], these data suggest the contribution 
of the p53 protein and associated pathways to the gen- 
eral cellular antiviral response against influenza infection 
and support the hypothesis that the virally induced glo- 
bal down-regulation of the p53 pathway leads to a cellu- 
lar state which favours their replication. 

Discussion 

A pathway of interest to study influenza interactions with 
the host cell 

In this present study, we report a microarray analysis 
revealing that infection by four human or avian influ- 
enza viruses (H3N2, H5N1, H5N2 and H7N1) signifi- 
cantly alters the gene expression of several host factors 



□ HCT116 p53-/- 
■ HCT116 p53+/+ 




24hpi 48hpi 

Figure 6 Comparison of viral titres (copies M/mL) produced in 
HCT116 p53 +/+ or -/- infected with H3N2 A/Moscow/1 0/99, 
calculated by RT-qPCR of M genome copies released in the 
supernatants. 



belonging to the p53 pathway. This pathway is com- 
posed of numerous genes, the protein products of which 
have a direct impact on host cell homeostasis, via the 
regulation of cell death and cell cycle progression as 
well as several other processes such as metabolism or 
the cellular immune response [23]. The transcriptional 
activity of p53 is tightly regulated by Hdm2 in response 
to several intracellular signalling cascades triggered by 
different stimuli (DNA damage, UV irradiation and 
hypoxia). The central core of the pathway includes the 
p53 protein itself, which can directly or indirectly shut 
down the growth of stressed cells mainly by inducing 
cell-cycle arrest and apoptosis. Considering the central 
role of the p53 pathway, the impact of influenza infec- 
tion on this cellular pathway is of particular interest. 

Several studies have underlined the interplay between 
influenza viruses and different cellular signalling path- 
ways, such as the PI3K/Akt, innate immune or apoptotic 
signalling pathways [4,36-39]. Intriguingly however, only 
a few have been dedicated to changes in p53 upon influ- 
enza infection, highlighting the need to further charac- 
terize the interplay between influenza viruses and the 
p53 pathway [25,26,34,40]. 

Influenza infection is associated with a global down- 
regulation of signalling cascades both upstream and 
downstream of p53 protein 

We first observed that upstream signals of the p53 path- 
way that directly connect stress (UV, hypoxia, DNA 
damage) -induced signalling cascades to regulators of 
p53, were generally down-regulated at the transcrip- 
tional level by infection. Among such regulating kinases, 
DNA-PK, Chk2, JNK1 and Gsk3/3 genes are all under- 
expressed following infection. 

Similarly, our results also revealed a significant 
decrease in gene expression for several endogenous p53 
transcriptional targets, including p21, 14-3-3, PERP, 
FAS, DR4/S, PIG3, BAX, Bcl-XL, PAI-1 and PCNA dur- 
ing infection by our set of viruses. Moreover, these 
results were confirmed at the protein level for p21, Bax 
and Bcl-XL. 

Altogether, our results indicate a global down-regula- 
tion in the mRNA expression of principal factors 
belonging to both upstream and downstream parts of 
the p53 pathway, suggesting a global negative effect on 
p53 pathway activity. 

Discrepancy between p53 mRNA, protein levels and 
activity in influenza infected cells 

The down-regulation of p53-target genes suggests that 
p53 expression and consequently its transcriptional 
activity is decreased during infection. We verified 
whether this was the case using luciferase assay and the 
monitoring of p53 phosphorylation status by western 
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blot in H3N2-infected A549 cells (data not shown) and 
revealed some discrepancies. On the one hand, no sig- 
nificant change in p53 mRNA expression was detected 
in infected cells, except for a significant decrease with 
H5N1. On the other hand, levels of p53 protein expres- 
sion were significantly increased by all the viruses, see- 
mingly incoherent with the reduced activity of p53 in 
infected cells. 

Altogether, these results suggest a likely complex reg- 
ulation of TP53 expression during infection. Regulating 
factors include transcriptional activators (ISGF3 PRKC8, 
HOXA5, CREBPB), inhibitors (BCL6, RBPJ) and stabiliz- 
ing factors (ELAVL1 and ZMAT3) [33]. Influenza infec- 
tion might therefore induce the over-expression of both 
inhibitors (BLC6 and RBPJ) and activators (ISGF3) (data 
not shown). In accordance with this notion, a recent 
study reported the stress-induced regulation of p53 
activity through the control of p53 mRNA stabilization 
and translation in addition to the well-described altera- 
tion of protein stability [33]. Other stabilizing/ destabiliz- 
ing factors, such as host microRNAs may also be 
involved and should be investigated in the context of 
infection [33], in particular with H5N1, which induces a 
down-regulation of p53 mRNA expression unlike the 
other viruses. 

Our results also clearly reveal that the impact that 
influenza viruses have on p53 is not exclusively on gene 
expression. Further investigations are necessary, in parti- 
cular on the stabilization and activation of p53 protein, 
considering the putative interaction of p53 with viral 
proteins during the time-course of infection. 

In this perspective, a recent study reported that the 
viral protein NS1 interacts with p53 and inhibits its 
activity [27]. Furthermore, transient expression of H5N1 
NS1 (wild-type or with mutations in regions 80-84 and 
92) reduced the transcriptional activity of p53 [41]. We 
can hypothesize that p53 protein is stabilized upon 
infection in an inactive form through interaction with 
the viral NS1 protein [27]. Such a scenario would 
explain our seemingly contradicting results concerning 
p53 expression and activity. This may even reflect what 
happens with a number of DNA viruses for which viral 
proteins like BZLF1 (Epstein Barr Virus) or HBx (Hepa- 
titis B virus) inhibit the transcriptional activity of a sta- 
bilized form of p53 [42,43]. Interestingly, the disruption 
of p53 signalling to p21 was described in human lung 
A549 cells in the case of BZLF1 viral protein [44]. 

Viral down-regulation of the p53 pathway and 
modulation of the cell cycle and apoptosis by influenza A 
infection 

We report the down-regulation of the p53-target genes 
p21 (Gl/S-arrest) and 14-3-3 (G2/M-arrest by sequester- 
ing Cdc2-complex Cyclin B), both of which are involved 



in cell cycle control. As expected, the down-regulation 
of p21 and 14-3-3 in association with an over-expression 
of the cell-cycle regulators TP53INP1 and GADD45G 
[35,45], could modulate cell cycle progression and 
resulted in Gl-arrest. Our results are in line with one 
recent study describing the effect of influenza A infec- 
tion on the cell cycle and show that viruses can induce 
an arrest in G0/G1 [46]. They also suggest a putative 
involvement of p21, TP53INP1 and Gadd45gamma in 
the influenza-induced cell cycle arrest. A thorough func- 
tional evaluation of these results is now needed which 
should contribute to a better understanding of the viral 
mechanisms involved which permit an optimal viral pro- 
tein expression and progeny production. 

The other group of down-regulated p53-target genes 
that we identified during infection encodes the pro- 
apoptotic proteins PERP [47], Fas, DR4/5, PIG3 and 
BAX. In addition, the Bax-inhibiting[48] anti-apoptotic 
factor Bcl-XL, was also under-expressed during infec- 
tion, as was PAI-1 whose role in regulating apoptosis 
remains unclear [49]. At the protein level, Bcl-XL 
expression was decreased by most viral infections. It is 
worth noting that the down-regulation of the three p53- 
target genes, p21, BAX and Bcl-XL, was considerably 
greater after H5N1 infection than with the other influ- 
enza A viruses tested. 

Concordantly, the knock-down of BAX was recently 
shown to reduce the replication of influenza virus [50]. 
Thus the down-regulation of a negative regulator of 
Bax, such as the anti-apoptotic Bcl-XL, at both mRNA 
and protein levels, such as we have reported here, would 
be expected to increase viral replication. Further silen- 
cing experiments are required to confirm such a result. 

Altogether, our observations suggest that influenza 
viruses may at least partially control intrinsic and extrin- 
sic pathways of apoptosis by decreasing the expression 
of both pro- and anti-apoptotic factors which are under 
the control of the p53 transcription factor, as suggested 
by several studies [39]. 

Is inhibition of the p53 pathway essential for the 
replication cycle of influenza viruses? 

In agreement with the above hypothesis, we showed 
that the level of H3N2 viral production is significantly 
higher in HCT116 p53 -/- than in HCT116 p53 +/+ 
cells, at both 24 and 48 hpi (figure 6). These results 
are in accordance with initial observations reported by 
Turpin and colleagues in a p53 dominant-negative 
A549 model and more recent data reported on siRNA 
high-throughput screening for host functional interac- 
tants [25,6]. 

These and our results suggest a marked antiviral facet 
for p53 and its pathway in the context of influenza 
infection, as has already been observed for several other 
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viruses. The antiviral effect of p53 depends largely on its 
pro-apoptotic activity, which limits virus replication 
[51]. However, other works have also demonstrated that 
p53 contributes to the innate cellular antiviral response 
by enhancing type I interferon (IFN)-dependent antiviral 
activity, independently of its function as a pro-apoptotic 
gene [52]. In light of our results on p53 targets impli- 
cated in cell-cycle regulation and a recent study on 
human circulating leukocytes from infected patients 
[53], the p53 antiviral role might involve not only its 
activity in mediating the IFN response and apoptosis, 
but also in controlling cell cycle progression. 

Conclusions 

In conclusion, the p53 pathway is inhibited during infec- 
tion by influenza A, either by the control over the 
expression of key p53 stabilizing factors, as is the case 
with H3N2, H5N1, H5N2 and H7N1, or by the decrease 
in p53 mRNA levels as is the case with H5N1. In accor- 
dance with an overall decrease of p53 activity, the global 
down-regulation of p53-target genes could contribute to 
changes in both the apoptotic response and cell-cycle 
progression, altogether inducing a pro-viral cellular con- 
text. The present study has highlighted the p53 pathway 
as an important player in influenza infection. Further 
investigations are now needed to discriminate the cellu- 
lar antiviral response from the viral hijacking of the p53 
pathway during the time-course of infection. 
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